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The two-body dynamics in general relativity has been solved perturbatively using the post- 
Newtonian (PN) approximation. The evolution of the orbital phase and the emitted gravitational 
radiation are now known to a rather high order up to 0(v 8 ), v being the characteristic velocity of 
the binary. The orbital evolution, however, cannot be specified uniquely due to the inherent freedom 
in the choice of parameter used in the PN expansion as well as the method pursued in solving the 
relevant differential equations. The goal of this paper is to determine the (dis)agreement between 
different PN waveform families in the context of initial and advanced gravitational-wave detectors. 
The waveforms employed in our analysis are those that are currently used by Initial LIGO/Virgo, 
that is the time-domain PN models TaylorTl, TaylorT2, TaylorT3, the Fourier-domain represen- 
tation TaylorF2 (or stationary phase approximant, SPA) and the effective-one-body (EOB) model, 
and two more recent models, TaylorT4 and TaylorEt. For these models we examine their overlaps 
with one another for a number of different binaries at 2PN, 3PN and 3.5PN orders to quantify their 
differences. We then study the overlaps of these families with the prototype effective-one-body fam- 
ily, currently used by Initial LIGO, calibrated to numerical relativity simulations to help us decide 
whether there exist preferred families, in terms of detectability and computational cost, that are the 
most appropriate as search templates. We conclude that as long as the total mass remains less than 
a certain upper limit M cr it, all template families at 3.5PN order (except TaylorT3 and TaylorEt) 
are equally good for the purpose of detection. The value of M cr it is found to be ~ 12 Mq for Initial, 
Enhanced and Advanced LIGO. From a purely computational point of view we recommend that 
3.5PN TaylorF2 be used below M cr it and EOB calibrated to numerical relativity simulations be used 
for total binary mass M > M cr it- 

PACS numbers: 



I. INTRODUCTION 



Sensitivity of several interferometric gravitational- wave detectors has either already reached, or is close to, the design 
goals that were set more than a decade ago [l], [H, H, 0, H, @, 01 ■ Upgrades that are currently underway and planned 
for the next four to five years, will see their sensitivity improve by factors of a few to an ordcr-of-magnitude 0]. 
Coalescing binaries consisting of neutron stars and/or black holes are probably the most promising sources for a 
first direct detection of gravitational waves. At current sensitivities, initial interferometers are capable of detecting 
binary neutron star inspirals at distances up to 30 Mpc, the range increasing to 60 Mpc for enhanced detectors (circa 
2009-2011) and 175 Mpc for advanced detectors (circa 2014+). Binary black holes or a mixed system consisting of 
a neutron star and a black hole can be detected to a far greater distance depending on the total mass and the mass 
ratio. 

The range of interferometric detectors for coalescing binaries is computed by assuming that one can pull the signal 
out of noise by matched filtering. This in turn means that one is able to follow the phasing of gravitational waves 
typically to within a fraction of a cycle over the duration of the signal in band. The reason for this optimism comes 
from the fact that one knows the phase evolution of the signal to a high order in post-Newtonian (PN) formalism 0] . 
Several authors have assessed whether the accuracy with which the formalism provides the waveforms is good enough 
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for the purpose of detection and parameter estimation [lfj, [ll|, QJ, |l3j, UA UM, UM LUX UM UM \M UM, \M, IM |25j . 
The problem, as we shall see below, is complicated since the PN approximation does not lead to a unique model 
of the phase evolution. Moreover, though PN results are good up to mildly rclativistic velocities, the standard PN 
approximants become less and less accurate in the strongly relativistic regime as one approaches the last stable orbit 
(LSO). Resummation methods [l4[ and in particular the effective-one-body (EOB) (2a. |27|. |28| extensions of the PN 
approximants, are needed for analytical treatments close to and beyond the LSO. 

The success in numerical-relativity simulations of binary black holes [2^, [3(| HH, H3, [H[ now provides results 
for gravitational waveforms that can be compared to standard PN results and other resummed extensions. On the 
one hand, the analytical PN results for the inspiral phase of the evolution are needed to calibrate and interpret the 
numerical relativity waveforms of coalescence and merger. On the other hand, the numerical relativity results extend 
the analytical approximations beyond the inspiral phase and provide the important coalescence and merger phases, 
producing the strongest signals that are crucial for the detection of binary black holes. However, numerical simulations 
are still computationally expensive and time-consuming and presently only a small region of the parameter space can 
be explored. Even in the foreseeable future, numerical relativity may not be able to handle, tens of thousands of 
cycles that are expected from highly asymmetric systems (e.g., a neutron star falling into an intcrmcdiatc-mass black 
hole of 100 Mq) or low-mass symmetric systems (e.g., a binary neutron star). Analytical models that smoothly go 
from the inspiral through coalcscense to quasi-normal ringing would be needed and this has led to phenomcnological 
templates (34 . l35l l36i . and EOB waveforms [HI, H3, [H, [H, [40, H3, E3, 0, EH ■ In particular, the recent improved 
EOB models [44], |45| which also incorporate a multiplicative decomposition of the multipolar waveform into several 
physically motivated factors supplemented by a suitable hybridisation (using test particle results) (HI, and an 
improved treatment of non-quasi-circular corrections, show evidence of remarkable success in modeling accurately the 
numerical relativity waveforms for different mass ratios. 

The emphasis of this work is different. Recently, there have been investigations [47j on the ability of various 



standard families of PN templates to detect a specific signal model TaylorEt [48|, I49I.I5C | and the often-used TaylorF2 
to detect a complete numerical relativity signal including merger and ringdown |35l. 36j. Reference [47[ modelled 
the signal by the TaylorEt approximant at 3.5PN order and looked at the cffcctualness and systematic biases in the 
estimation of mass parameters for TaylorTl, TaylorT4 and TaylorF2 templates in the LIGO and Virgo detectors. It 
also looked into the possibility of improving the effectualness by using unphysical values of v beyond the maximum 
value of 0.25. It was found that the overlaps of a TaylorEt signal with TaylorTl, TaylorT4 and TaylorF2 template 
is smaller than 0.97 and involved for equal-mass systems a large bias in the total mass. For unequal-mass systems 
higher overlaps can be obtained at the cost of a large bias in mass and symmetric mass ratio v and which can be 
further improved by unphysical values of v > 0.25. The templates arc more unfaithful with increasing total mass. 
To detect optimally the complete numerical-relativity signal, including merger and ringdown, Ref. [35J suggested the 
possibility of using the TaylorF2 template bank with a frequency cutoff f c larger than the usual upper cutoff (i.e., the 
Schwarzschild LSO) and closer to the fundamental quasi-normal mode frequency of the final black hole. Moreover, 
they proposed to further improve this family by allowing either for unphysical values of v or for the inclusion of a 
pseudo 4PN (p4PN) coefficient in the template phase, calibrated to the numerical simulations. Reference [36[ extended 
the results of Ref. [35l | to more accurate numerical waveforms, found that 3.5PN templates are nearly always better 
and rarely significantly worse than the 2PN templates, and proposed simple analytical frequency cutoffs for both 
Initial and Advanced LIGO for example for Initial LIGO they recommended a strategy using p4PN templates for 
M < 35M Q and 3.5PN templates with unphysical values of v for larger masses. However, we notice that there is no 
reason for changing the template bank above 35M . Reference [3|| could have used the p4PN templates over the 
entire mass region, if they had not employed in their analysis the p4PN coefficient used in Ref. (35[ , but had calibrated 
it to the highly accurate waveforms used in their paper 1 . 

In this work our primary focus is on binary systems dominated by early inspiral and on a critical study of the 
variety of approximants that describe this. Towards this end, in this paper we will provide a sufficiently exhaustive 
comparison of different PN models of adiabatic inspiral for an illustrative variety of different systems and quantify how 
(dis)similar they are for the purpose of detection. The choice of the PN models used in this paper is motivated by the 
fact that they are available in the LIGO Algorithms Library (LAL) and some of them have been used in the searches 
by Initial LIGO. We also compare all these PN models with one fiducial EOB model calibrated to numerical- relativity 
simulations (39j to delineate the range of mass values where one must definitely go beyond the inspiral-dominated 
PN models to a more complete description including plunge and coalescence. The choice of this fiducial, preliminary 
EOB model is only motivated by the fact that it is the EOB model available in LAL and it is currently used for 
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searches by Initial LIGO. It will be improved in the future using the recent results in Refs. [HI, iH|. We will conclude 
that for total masses below a certain upper limit M cr it; all template families at 3.5PN order (except for TaylorT3 
and TaylorEt) are equally good for the purpose of detection. M cr it is found to be ~ 12 Mq for Initial, Enhanced and 
Advanced LIGO. Based solely on computational costs, we recommend that 3.5PN TaylorF2 be used below M cr it and 
EOB calibrated to numerical relativity simulations be used for total binary mass M > M cr [ t . 

The paper is organized as follows. In Sec. [ill we summarise the present status of the PN approximation. In Sec. IIIII 
we recapitulate for completeness the main PN approximants and try to provide a ready-reckoner for the equations 
describing them and the relevant initial and termination conditions. In Sec. IIVI we discuss the frequency evolution 
in each of these models. In Sec. |V]we discuss overlaps and the maximization used in this work. Section [VTI and IVTll 
presents the results of our analysis related to the effectualness, while Sec. IVIIII summarizes the results related to the 
faithfulness. In Sec. ILXI we summarize our main conclusions. Readers who are interested in the main results of the 
paper and want to avoid technical details could skip Sees. HH IIIII IIVI and[V] read the main results of Sees. IVII IVIII 
and IVIlfl and mainly focus on Sec. ILXI 



II. CURRENT STATUS OF POST-NEWTONIAN APPROXIMATION 



Post-Newtonian approximation computes the evolution of the orbital phase 4>{t) of a compact binary as a pertur- 
bative expansion in a small parameter, typically taken as v — (irMF) 1 / 3 (characteristic velocity in the binary), or 
x = v , although other variants exist. Here M is the total mass of the binary and F the gravitational- wave frequency. 
In the adiabatic approximation, and for the restricted waveform in which case the gravitational wave phase is twice 
the orbital phase, the theory allows the phasing to be specified by a pair of differential equations <fi(t) = v 3 /M, 
v = — T (v) / E' (v) , where M is the total mass of the system, T its gravitational-wave luminosity and E'(v) is the 
derivative of the binding energy with respect to v. Different PN families arise because one can choose to treat the ratio 



J-/E'(v) differently while being equivalent with the same PN order 17[. For instance, one can leave the PN expansions 
of the luminosity T(v) and E'(v) as they appear (the so-called TaylorTl model), or expand the rational polynomial 
T(v)/E'{v) in v to consistent PN order (the TaylorT4 model), recast as a pair of parametric equations 4>(v) and t(v) 
(the TaylorT2 model), or the phasing could be written as an explicit function of time cf>(t) (the TaylorT3 model). 
These different representations are made possible because one is dealing with a perturbative series. Therefore, one is 
at liberty to "resum" or "reexpand" the series in any way one wishes (as long as one keeps terms to the correct order 
in the perturbation expansions), or even retain the expression as the quotient of two polynomials and treat them 
numerically. There is also the freedom of writing the scries in a different variable, say (suitably adimensional) E (the 
so called TaylorEt model). 

In addition to these models, there have been efforts to extend the evolution of a binary beyond what is naturally 
prescribed by the PN formalism. Let us briefly discuss two reasons why the PN evolution cannot be used all the way 
up to the merger of the two bodies. PN evolution is based on the so-called adiabatic approximation according to which 
the fractional change in the orbital frequency F or b over each orbital period is negligibly small, i.e. F or \,/F^ Tb <C 1. 
This assumption is valid during most of the evolution, but begins to fail as the system approaches the LSO where 
/lso = (6 3 / 2 7rM) -1 . In some cases, the frequency evolution stops from being monotonic and / changes from being 
positive to negative well before reaching the LSO — an indication of the breakdown of the approximation. 

From the view point of maximizing detection potential one is also interested in going beyond the inspiral phase. 
The merger and ringdown phases of the evolution, when the luminosity is greatest, cannot be modelled by standard 
PN approximation. The use of resummation techniques more than a decade ago was followed by the construction 
of the EOB model [H, [U HH, which has analytically provided the plunge, merger and ringdown phases of the 
binary evolution. As mentioned before, more recently, these models have been calibrated to numerical relativity 
simulations [H, H3, [H, [H, HO, SH, S3, EE 0, Hf| . We now have a very reliable EOB model that can be used to model 
the merger dynamics. 

An astronomical binary is characterized by a large number of parameters some of which are intrinsic to the system 
(e.g., the masses and spins of the component stars and the changing eccentricity of the orbit) and others that are 
extrinsic (e.g., source location and orientation relative to the detector). In this paper we will worry about only the 
detection problem. Furthermore, we will assume that a coincident detection strategy will be followed so that we do not 
have to worry about the angular parameters such as the direction to the source, wave's polarization, etc. If binaries 
start their lives when their separation r is far larger compared to their gravitational radius (i.e., r 3> GM/c ), by 
the time they enter the- sensitivity band of ground-based detectors any initial eccentricity would have been lost due 
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to gravitational radiation reaction, which tends to circularize 2 a binary [58|, (59|. Therefore, we shall consider only 
systems that are on a quasi-circular inspiralling orbit. We shall also neglect spins which means that we have to worry 
in reality about only the two masses of the component bodies. 

Our goal is to explore how (dis)similar the different waveform families are. We do this by computing the (normalized) 
cross-correlation between signals and templates, maximized either only over the extrinsic parameters of the templates 
(faithfulness) or over the intrinsic and extrinsic parameters of the templates (effectualness) , the noise power spectral 
density of the detector serving as a weighting factor in the computation of the correlation (see SeclVj). Our conclusions, 
therefore, will depend on the masses of the compact stars as well as the detector that we hope to observe the signal 
with. 

The overlaps (i.e., the normalized cross-correlation maximized over various parameters and weighted by the noise 
power spectral density) we shall compute are sensitive to the shape of the noise spectral density of a detector and 
not on how deep that sensitivity is. Now, the upgrade from initial to advanced interferometers will see improvements 
in sensitivity not only at a given frequency but over a larger band. Therefore, the agreement between different PN 
models will be sensitive to the noise spectral density that is used in the inner product. Thus, we will compare the PN 
families using power spectral densities of initial and advanced interfcrometric detectors. 

We end this brief overview with the following observation. As mentioned earlier, following all present gravitational 
wave data analysis pipelines, this paper works only in the restricted wave approximation. This approximation assumes 
the waveform amplitude to be Newtonian and thus includes only the leading second harmonic of the orbital phase. 
Higher PN order amplitude terms bring in harmonics of the orbital phase other than the dominant one at twice the 
orbital frequency. Their effects can be significant [6(J [HI, especially close to merger [44| . and they need to be 
carefully included in future work. 



III. THE PN APPROXIMANTS 



For the convenience of the reader, in this section, we recapitulate the basic formulas for the different PN families 
from Refs. [13, El- While comparing the expressions below to those in Refs. [13, El recall A = -1987/3080 [H[63j| 
and O = —11831/9240 [64|,[65j|. In addition to the evolution equations, we shall also provide initial and final conditions. 
From the perspective of a data analyst, the initial condition is simply a starting frequency Fq and phase <po, which 
can be translated, with the help of evolution equations, as conditions on the relevant variables. We shall also give 
explicit expressions for the evolution of the gravitational wave frequency, namely F = dF/dt, or more precisely, the 
dimensionless quantity F F~ 2 , in Sec. IIV1 where they will be used to study the rate at which the binary coalesces 
in different PN families, which will help us understand the qualitative difference between them. The contents of this 
section should act as a single point of resource for anyone who is interested in implementing the waveforms for the 
purpose of data analysis and other applications. 

The basic inputs for all families are the PN expressions for the conserved 3PN energy (per unit total mass) [62|, 
[H, [H, [|3, H M E 3(v) and 3.5PN energy flux [H, M, M, IzE Izl ^3.s(v), 
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where 7 = 0.577216 ... is the Euler constant. In the adiabatic approximation one assumes that the orbit evolves slowly 
so that the fractional change in the orbital velocity ui over an orbital period is negligibly small. That is, — <C 1, 



2 Though this assumption is justified for the prototypical binaries we focus on in this work, there exist credible astrophysical scenarios that 
lead to inspiral signals from binaries with non-ne glig ible eccentricity in the sensitive detector bandwidth. A more involved treatment is 
then called for and available. See e.g. [H [H, |M 111, HI, HI El • 
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or, equivalcntly, -jjj <C 1. In this approximation, one expects the luminosity in gravitational waves to come from 
the change in orbital energy averaged over a period. For circular orbits this means one can use the energy balance 
equation T = —dE/dt where £ = ME. 

In the adiabatic approximation one can write an equation for the evolution of any of the binary parameters. For 
instance, the evolution of the orbital separation r(t) can be written as r(t) = £/{d£/dr) = —J-/{d£/dr). Together 
with the Kepler's law, the energy balance equation can be used to obtain the evolution of the orbital phase 3 : 



dt M 

dv 

~dt + ME'(v) 

or, equivalcntly, 



= 0, (3.3a) 
= 0, (3.3b) 



t(v) = tref + M / dv—j-^-, (3.4a) 



frof 



,3 E'(v) 



(v) = 4> Te{ + / dvv A —^-, (3.4b) 



where i re f and </> re f are integration constants and v le { is an arbitrary reference velocity. 



A. TaylorTl 

The TaylorTl approximant refers to the choice corresponding to leaving the PN expansions of the luminosity J-(v) 
and E'(v) as they appear in Eq. (|3.3p as a ratio of polynomials and solving the differential equations numerically 



(Tl) 



dt M 

dv 

~dt + ME'(v) 



= 0, (3.5a) 
= 0. (3.5b) 



In the above v = z/ T1 ) but for the sake of notational simplicity we write only v; from the context the meaning should 
be clear. In the formulas of this section, and in the sections that follow, the ex pre ssions for J-(v) \E(v)\ are to be 
truncated at relative PN orders 2 [2], 3 [3] and 3.5[3] to obtain 2PN [U [z!, H El , 3PN and 3.5PN [H, [13, [z3] 
template or signal models respectively. 

To see how to set up initial conditions, refer to Eq. (|3.4[) . Let the initial gravitational wave frequency be Fq or, 
equivalently, initial velocity vo = {-kMFq) 1 ^ . One normally chooses t = at v = Vq- This can be achieved by choosing 
^rcf = vq and i re f = 0, in Eq. (|3 .4|) . The initial phase (f> re { is chosen to be either or tt/2 in order to construct two 
orthogonal templates (see Sec. IV Al for details). 



B. TaylorT4 

Taylor T4 was proposed in Ref. [22j and investigated in Refs. [13, H3, , thus many years after the other approxi- 
mants discussed in this paper were proposed (with the exception of TaylorEt, which is even more recent). However, it 
is a straightforward extension of TaylorTl and at 3.5PN order by coincidence is found to be in better agreement with 
numerical simulations of the inspiral phase [32l . [35l . |3?1 140| . l42l . |48| . |76| . The approximant is obtained by expanding 
the ratio of the polynomials J-{v)/E'(v) to the consistent PN order. The equation for t/ T4 )(i) = v{t) at 3.5PN order 
reads, 



3 Recall that the gravitational-wave phase is twice the orbital phase for the restricted waveform and leads to differences in factors of 2 
between the equations here for the orbital phase and those in llTt for the gravitational-wave phase. 
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The orbital phase (p^ 14 ^ is determined, as in the case of TaylorTl, by Eq. (|3.3ap and numerical solution of Eq. (|3.6p 
and (|3.3a|) yields the TaylorT4 approximant. 

Note that although TaylorTl and TaylorT4 are perturbativcly equivalent, the evolution of the phase can be quite 
different in these two approximations. The asymptotic structure of the approximants are also quite different: while v 
can have a pole (although not necessarily in the region of interest) when using Eq. (|3.5b[) none is possible when Eq. 
(|3.6[) is used. Differences of this kind can, in principle, mean that the various PN families give different phasing of 
the orbit. The hope is that when the PN order up to which the approximation is known is large, then the difference 
between the various PN families becomes negligible. 

Setting up the initial conditions for TaylorT4 is the same as in the case of TaylorTl. 



C. TaylorT2 

TaylorT2 is based on the second form of the phasing relations Eq. p.4[) . Expanding the ratio of the polynomials 
T(v)/E'{v) in these equations to consistent PN order and integrating them one obtains a pair of parametric equations 
for <j)(v) and t(v), the TaylorT2 model. 
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Of all models considered in this study, TaylorT2 is computationally the most expensive. This is because the phase 
evolution involves solving a pair of transcendental equations which is very time-consuming. 
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In this case, t Ie f has to be chosen so that t = when F = F or v = Vq. This can be achieved most simply by 
solving for i re f, using Eq. (|3.8b[) . substituting v — vq on the right hand side and putting the left side to zero. 
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D. TaylorT3 

This form of the approximant goes a step further than the previous TaylorT2 approximant. After computing as 
before a parametric representation of the phasing formula <p(v) and t(v), one explicitly inverts t(v) to obtain v(t) and 
uses it to produce an explicit representation of 4>{t) = <fi(v(t))). This is the TaylorT3 approximant: 

n 

€/2(t) = ^S 3) +^E^ fc ' ( 3 - 9a ) 

fc=0 

n 

^/2 3) W = FnT,^ 6 ^ ( 3 - 9b ) 

k=0 

where 6 = [^(tref — i)/(5M)] -1 / 8 and F = (2d(j)/dt){2^)^ 1 = v 3 /(ttM) is the instantaneous gravitational-wave 
frequency. 
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The initial conditions in this case is slightly more complicated than the previous cases. Given an initial frequency 
Fq, one numerically solves Eq. (|3.10b[) to find the value of t Ie f at which F — Fq and t = (recall that involves t re f.) 
Note that as t — > t le {, formally F — > diverges. 



E. TaylorEt 

The TaylorEt was recently introduced in Ref. [48|, [4^, [5C| ■ Introducing 4 £ = —2E/v (recall that our E is conserved 
energy per total mass), the TaylorEt approximants are obtained starting from Eq. (|3.ip for E(x) or (,{x) and inverting 
it to obtain x(Q): 



x = C 



1 
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(3.11) 



With this choice of variable the equation determining the evolution of v, Eq. (|3.3b[) , transforms to the balance equation 
for E rewritten in terms of the C variable: 



d( = 2T(v(Q) 
dt v M 



(3.12) 



There is no difference between Tl and T4 approximants in the Et-parametrisation and the gravitational-wave phasing 
equations Eq. (|3.3a|) and Eq. (|3.3bj) in terms of £ become [13] , 



4 Note that the £ in this paper is denoted variously by £ in [49ll but by £ in e.g. [4?]]. 
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(3.13b) 



To set up the initial condition note that 2ttF = 2 d(f>/dt. Given an initial frequency F one finds the initial value £ of 
C by numerically solving Eq. (|3.13a|) . by setting the left hand side to ttFq. 



F. TaylorF2 

The most commonly used form of the approximant is the Fourier representation computed using the stationary 
phase approximation (SPA). Using the SPA the waveform in the frequency domain may be written as, 



/* spa (/) 



Htf) 



4, f {t) = 2nft-24>{t), 



(3.14) 



where tf is the saddle point defined by solving for t, dipf(t)/dt = 0, i.e. the time tf when the gravitational-wave 
frequency F(t) becomes equal to the Fourier variable /. In the adiabatic approximation, (denoting vj = (ttM/) 1//3 ) 
the value of tf and ipf(tf) are given by the following integrals: 



trei + M 

2tt ft le { - 



dv, 



cf + 2 



dv. 



As in the time domain case it is more efficient to use the equivalent differential form 

ttM 2 E'{f) 



dib dt 

i- 2 ^^ d f 



3« 2 F(f) 



0, 



(3.15a) 



(3.15b) 



(3.16) 



and this characterizes the TaylorFl approximant. 

The analogue of the TaylorT2 in the frequency domain follows by explicitly truncating the energy and flux functions 
to consistent post-Newtonian orders and explicating the v- integration in the above. This leads us to a Fourier domain 
waveform, the TaylorF2, which is the most often employed PN-approximant, given by 

Hf) =AT 7/6 e W) , (3.17) 

where A oc A^ 5 / 6 Q(angles)/£>, and D the distance to the binary. To 3.5PN order the phase of the Fourier domain 
waveform is given by 



^ F 5 2) (/) 
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where v = (irMf) 1 / 3 . 

In this case one has to specify the constants t c and <f> c and they can be chosen arbitrarily. 



G. The effective-one-body model 



In this paper since we are not particularly concerned with the coalescence signal, we employ the less sophisticated 
earlier version of the EOB model calibrated to numerical-relativity simulations from Ref. [39| (f° r more sophisticated 
versions of the EOB model see Refs. 0, HH, |H, H3, IH). Below we briefly review the EOB model from Ref. [39] . 

Introducing polar coordinates (r, <fi) and their conjugate momenta (p r ,P(j>), the EOB effective metric takes the 
form [H 



dsi 



-A(r) dt 2 



D(r) 
Mr) 



dr A 



d6 z 



The EOB Hamiltonian reads 



H ieaX {r,p r ,pt) = ^H rcal = M J I + 2v 
with the effective Hamiltonian [20, u8f 



H oit - n 



H eS (r, Pr ,pt) = (,H eS = ix 



\ 



Air) 



1 



Mr) 
D(r) 



Pl + -i+ 2(4 - 3v) v 



Pi 



The Taylor-approximants to the coefficients A(r) and D(r) can be written as [2a. \2i 



fe+i 



Mr) - Y.^ 

i=0 



(3.19) 
(3.20) 

(3.21) 

(3.22a) 
(3.22b) 



The functions A(r), D(r), A^{r) and D^{r) all depend on the symmetric mass ratio v through the ^-dependent 
coefficients aj(v) and di(v). These coefficients arc currently known through 3PN order (i.e. up to k = 4) and can be 
read from Ref. [39|. During the last stages of inspiral and plunge 5 , the EOB dynamics can be adjusted closer to the 
numerical simulations by including in the radial potential A(r) a p4PN coefficient 05 (y) and a^(i>) = Xq u, with A a 
constant 6 . In order to assure the presence of a horizon in the effective metric (|3.19j) . a zero needs to be factored out 
from A(r). This is obtained by applying a Pade resummation (28|. The Pade coefficients for the expansion of A(r) 
and D(r) at p4PN order are denoted A\(r) and D®(r), and their explicit form can be read from Ref. (39j . 

The EOB Hamilton equations are written in terms of the reduced (i.e., dimensionlcss) quantities jj real [defined in 
Eq. (ESP ], t = t/M, and Q = uM |27j]: 



dr 
dt 
d0 
dt 

dp r 

dt 
dp± 
dt 



dH 



real 



dp r 

8H real 



{r,Pr,P4>) , 
(r,Pr,P4,) , 



dH 



real 



dr 

F<t>{r,Pr,P<t>) , 



(r,p r ,p<t>) , 



(3.23a) 

(3.23b) 

(3.23c) 
(3.23d) 



5 To deal with the steep rise of various quantities during the plunge, it is advantageous to consider the EOB equations in terms of the 
tortoise radial coordinate r„ and its conjugate p r „ rather than in terms of the standard radial coordinate r and p r as above. The form 
of H eB in the two cases will be different [38j [. For the level of accuracy in our present work, this difference is irrelevant. 

6 Note that A was denoted A in Ref. [H, and a 5 in Refs. [38l , liM . 113 . itl . 
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with the definition u) = dxj>/dt. Another critical input to the EOB model is the form for the radiation reaction 
force arising from the basic PN expression of the energy flux. Different choices include Pade resummation [Til ], and 
the more recent pi m - resummation |46| . It also further includes the introduction of terms describing next-to-quasi- 
circular effects. Here, for the component of the radiation-reaction force we use the less sophisticated Keplerian 
Pade-approximant to the energy flux as given by Eq. (15) of Ref. [39l | . 

The inspiral-plunge EOB waveform at leading order in a PN expansion reads 

^insp-plunge^ = -1/3 cos [ 2 0(i)] . (3.24) 

The merger-ringdown waveform in the EOB approach is built as a superposition of quasi-normal modes [27|, H3, [H, 
H, SO, [231 , as 

N-l 

h mersei - RB (t) = A n e^** - *™"* 5 , (3.25) 

n=0 

where n is the overtone number of the Kerr quasi-normal mode, N is the number of overtones included in our 
model, and A n are complex amplitudes to be determined by a matching procedure described below. The quantity 
er„ = cu n — ia n , where the oscillation frequencies ui n > and the inverse decay-times a n > 0, are numbers associated 
with each quasi-normal mode. The complex frequencies are known functions of the final black-hole mass and spin 
and can be found in Ref. (78j . The final black-hole masses and spins are obtained from the fitting to numerical results 
worked out in Ref. 39] . 

The complex amplitudes A„ in Eq. (|3.25j) are determined by matching the EOB merger-ringdown waveform with the 
EOB inspiral-plunge waveform close to the EOB light ring. In particular, in Ref. (39j the matching point is provided 
analytically by Eq. (37). In order to do this, N independent complex equations are needed. The N equations are 
obtained at the matching time by imposing continuity of the waveform and its time derivatives, 

^in S p-plu„gc (tmatch) = ^mcrgcr-RD ( £ match) ^ (fc = 0, 1, 2, • • • , N — 1) . (3.26) 

In this paper we use N=3. The above matching approach is referred to as point matching. It gives better smoothness 
around the matching time, but it is not very stable numerically when N is large and higher order numerical derivatives 
are needed. More sophisticated matching procedures have been proposed in the literature to overcome the stability 
issue. Reference (38| introduced the comb matching approach where N equations are obtained at N points evenly 
sampled in a small time interval Ai matc h centered at t ma tch- More recently, to improve the smoothness of the comb 
matching Ref. [45| introduced the hybrid comb matching where one chooses a time interval Ai matc h ending at i ma tch, 
and imposes not only the continuity of the waveform at N — 4 points evenly sampled from f roa tch — A£ matc h to i ma t c hj 
but also requires continuity of the first and second order time derivatives of the waveform at i ma t c h — Ai matc h and 

^match- 

Finally, the full (inspiral-plunge-merger-ringdown) EOB waveform reads 

hit) = ftto-P-plunge^) 0( tmatch _ t ) + ^merger-RD ^ _ ^ (3 2?) 

where we denote with 6 the Heaviside step function. 

H. Waveforms and termination conditions 

Before concluding this Section we note a few other points concerning the generation of the waveform. Since our 
goal is to study the agreement between different waveforms it is not necessary to separately consider the two different 
polarizations but only the detector response. For time-domain models TaylorTl, TaylorT2, TaylorT3, TaylorT4 and 
EOB the waveform is taken as: 

h A (t) = Cv% 8in[2 A (t)], 

where va and <f>A (i) are computed using the relevant formulas corresponding to the approximant A. In the case of 
TaylorEt the waveform is taken to be 

h Et (t)=C{(t) sin[2 0Et(i)]- 

In all cases the constant C is fixed by demanding that the norm of the signal be unity (cf. Sec.|V|. The initial phase 
of the signal is set to 0, while in the case of templates we construct two orthonormal waveforms corresponding to the 
starting phases of and tt/2. 
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TABLE I: Termination condition for waveform generation is chosen to be either LSO corresponding to Schwarzschild metric 
vs — 6 _1//2 , or the extremum defined by the P-approximant of the energy function as in [lj which is «p 4 at 2PN and Vp g at 3- 
and 3.5PN. In the case of TaylorT3 at 3.5PN, as the frequency evolution is not monotonic, the evolution has to be terminated 
prematurely at v m such that F(v m ) = 0. 



Order/Approx 


Tl 


T2 


T3 


T4 


Et 


F2 


2PN 


vs 


VS 


VS 


Vp 4 


vp 4 


vp 4 


3PN 


VS 


Vs 


Vs 


«P 6 


VP 6 


vp 6 


3.5PN 


VS 


vs 


Vm 


vp 6 


VP 6 


vp 6 



The waveforms are terminated when v reaches the value quoted in Table I or before, if the frequency evolution 
is not monotonic (see next Section). For instance, in the case of TaylorT3 at 3.5PN order the approximant has an 
unusual behaviour whereby the frequency evolution ceases to be monotonic well before v reaches the nominal value of 
1/V6. In the case of TaylorTl, TaylorT2 and TaylorT3, the termination is at the LSO defined by the Schwarzschild 
metric, namely v = l/y/E, at all PN orders, but we also check for monotonicity of the frequency evolution. For other 
approximants, except EOB, we terminate at the extremum of the P-approximant energy function [bH ]. In the case of 
EOB, the waveform is terminated at the end of the quasi-normal ringing. 



IV. FREQUENCY EVOLUTION 



The quantity that determines the evolution of a binary, its phasing and the duration for which it lasts starting from 
a particular frequency, is the acceleration of the bodies under radiation reaction. Equivalently, it is the evolution of 
the derivative of the gravitational wave frequency F = dF/dt, which determines the phasing of the waves. When 
the separation between the bodies is large, the frequency evolution is slow and the quantity [l~7| e(t) = FF~ 2 , which 
measures the fractional change in the frequency over a period, is small: FF~ 2 <C 1. As the binary evolves, this 
quantity increases but, as seen in numerical evolutions, remains finite and positive all the way up to the merger of 
the two bodies. In what follows we will explore the behaviour of e as a function of the PN parameter v rather than 
t, because the former parameter is (mass) scale free, unlike the latter. 

Computing the adiabaticity parameter e(v) in the case of TaylorTl and TaylorT4 is straightforward using Eqs. 
(|3.5b[) and (|3.6|) . In the case of TaylorT2, one differentiates Eq. (|3.8b[) with respect to v and then takes its reciprocal. 
Finding e(v) in the case of TaylorEt is more involved. The frequency F is given by Eq. (|3.13a[) but the right hand 
side is a function of £. One must, therefore, combine Eqs. (|3.13a[) and (|3.13b() to find the derivative of the frequency: 

dC, dt d( \dt J dt y ' 

The above equation still gives F as a function of £■ One can then use Eq. p. lip to get e(v). Consequently, there 
is no guarantee that v will be monotonic in the region of interest. However, we do find that the function eEt(^) is 
positive in the region of interest and therefore v increases monotonically for TaylorEt. To find e(v) for TaylorT3, F 
is given by differentiating Eq. (|3.10b[) with respect to t (recall = 6{t)) and then one uses the same equation to find 
v = (ttMF) 1 / 3 at a given t. It turns out that for TaylorT3 the function ex3 can become negative in the region of 
interest (exactly when this happens depends on the PN order and mass ratio) and so v does not generally increase 
monotonically. 

Figure[TJ left panel, plots e(v) for two values of the mass ratio: v = 0.10 and v = 0.25. When v is small (y <C 1/V6) 
e(v) for the different approximants is the same. Therefore, in systems for which v remains small when the signal is 
in band (as, for example, in a binary neutron star), the different approximants, as we shall see in the next Section, 
agree well with each other. As v approaches l/v6j different approximations tend to differ greatly, which means we 
cannot expect good agreement between the different PN families. Of the approximants considered here, TaylorEt 
seems to have the smallest value of e(v) at any given v. Therefore, the evolution will be slower, and the duration of 
the waveform from a given frequency larger, than the other approximants [5(J. TaylorT3 also differs from all others 
because e(v) becomes negative before the last stable orbit, and so v does not generally increase monotonically for this 
approximant. This behavior can be seen at 2PN and 3.5PN orders in the left panel of Fig.[T] The reason for this can 
be seen in Fig. [T] right panel, where we have shown the time development of £t3(£) for two values of v = 0.10, 0.25. 
Since F becomes negative before reaching the last stable orbit, the waveform has to be terminated before before v 
reaches l/-\/6. 
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FIG. 1: On the left hand panel the plots show the evolution of frequency in different PN families. The adiabaticity parameter 
e(t) = F~ 2 F is essentially the same for all the different approximations at « < 1 As the binary gets close to coalescence the 
various approximations begin to differ from each other. The right hand panel shows the adiabaticity parameter as a function 
of time t at 3.5PN order. Note that 6t3(£) begins to decrease and even becomes less than zero before v reaches its nominal 
value of l/y/6. This leads to waveforms that are significantly shorter in the case of TaylorT3. 



/ieob(xo) 




FIG. 2: Schematic plot of distance (or mismatch) relation between templates and exact, numerical and EOB waveforms. 



V. EFFECTUALNESS 



The goal of this study is to compare the different PN approximations by measuring their mutual effectualness (i.e., 
overlaps maximized over intrinsic and extrinsic parameters) for a number of different mass pairs. To this end it will 
be very useful to define the scalar product of waveforms. Given waveforms hk and q/-, k = 0, . . . , N — 1, where hk 
is the kth sample of the signal h(t) at time tk = fcA, A = l/f s being the sampling interval corresponding to the 
sampling rate f s , their scalar product is defined by 7 

N-l A 

(h, q) (r k ) = 2 £ [H m Q* m + ff* Q m ] e -«/iv (51) 



It is conventional to define the scalar product in the continuum limit. Here, however, we have given the definition for a discretely 
sampled data and this is the expression that is used in computing the overlaps. 
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where Af = f s /N, f m = mA/, Tk = kA is the lag of the template a measure of the relative time-shift between the 

template and signal, H m = A^^Jq 1 hk e 2lTlmk / N is the discrete Fourier transform of h(t) (similarly, Q m ) and Sh{fm) 
is the one-sided noise power spectral density of a detector. In comparing two waveforms the overall amplitude is of 
no interest and we should, therefore, consider waveforms with unit norm, namely h = hj yj (h, h). Consequently, the 
relevant quantity is the scalar product between normalized waveforms defined by 

0[h, q] = (h, q) (5.2) 



A. Maximization of the overlaps 



The signal and the template both depend on a set of parameters of the source (e.g. , masses and initial spins of the 
component masses) and its orientation relative to the detector. We shall be concerned with binaries with non-spinning 
components on quasi-circular orbits. Such systems are characterized by two intrinsic parameters, namely the masses 
mi and m 2 of the components, and two extrinsic parameters, namely the time-of-coalescence tc and the phase of the 
signal at that time 4> c . The overlap integral, therefore, depends on the parameters of the signal and the template and 
the relevant quantity is the overlap maximized over these parameters. 

The data analysis problem is concerned with digging out a specific signal buried in noisy data. This means that the 
parameters of the signal are fixed but the data analyst is at liberty to maximize over the parameters of the template. 
In this paper we will explore the cffcctualncss of templates; that is to say the overlap maximized over a template's 
parameters keeping those of the signal fixed. We will do this for several choices of the component masses of the binary. 
However, the time-of-coalescence tc and the phase 4>c of the signal at that time, are arbitrarily chosen to be equal 
to zero. A caveat is in order concerning the value of the effectualness arising as a result of our choice of tc and <j>c ■ 
the maximized overlap is not very sensitive to our choice of tc but it could vary by several percents depending on the 
choice of a signal's phase, especially when the signal and the template families are not very close in the geometrical 
sense. 

Maximization over a template's masses is carried out using a bank of templates and the template bank is set up 
such that for all signals of the same family as the template their best overlap with the nearest template is larger than 
a certain value called the minimal match MM. Our template placement is as in Ref. [79l |. which is known to produce, 
with probability close to 1 matches larger than the minimal match for the TaylorTl, TaylorT3, TaylorF2 and 
EOB families of signals (and templates) for the range of masses considered in this paper. We have checked this to be 
true also for TaylorEt and TaylorT4 families. 

We have used a minimal match of MM = 0.99 in all cases. Maximization over time-of-coalescence is accomplished 
by looking at the overlap integral at different lags Tk- Finally, since our templates are of the form hk = Ak cos(4>k + 4'o)t 
where 4>o is an unknown constant phase offset, maximization over tfio can be achieved by using two quadratures of the 
template, h° k = Ak cos(</>fc) and ti*J 2 = Ak cos(</>fc + 7r/2) : 



maxO[/i, q] = 

<Po 



(h°, qf + (/W 2 , g y 



1/2 

(5.3) 



When the signal and the template belong to the same family the maximized overlap is at least MM. When the 
waveforms belong to different families the maximized overlap is less than MM. 

Our approach to finding the effectualness of a template with a signal of "fixed" parameters is here somewhat different 
from what is normally followed in the literature, but more appropriate in the context of data analysis. In the literature 
on the comparisons of different PN models, one normally measures either the best or the minimax overlap [Til ]. The 
best overlap gives the maximum of the overlap over the masses and tc but maximized over the constant phases of 
both the signal and the template. On the other hand, the minimax overlap is the overlap maximized over the masses 
and tc but minimized over the constant phases of the signal and the template. As mentioned earlier, we fix the 
phase of the signal to be equal to zero and hence our effectualness is, in principle, smaller than best overlaps but 
larger than minimax overlaps. The difference between the best and minimax overlaps is tiny when the effectualness is 
intrinsically large (i.e., close to 1), but could differ by 5 — 8% when the best overlap is ~ 0.8. This should be kept in 
mind while interpreting our results. Moreover, as mentioned earlier, instead of numerically searching for the maxima 
of the overlap in the space of masses we just use a grid of templates with a minimal match of MM = 0.99. 

We will compute effectualness between every possible template and signal. If our template is the PN approximation 
A and the signal is the PN approximation B then we are interested in computing the matrix eab defined by 

e AB = m a xO[h A (\ A ), h B (X B )}, (5.4) 
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where X A and X B are the parameters of the template and the signal, respectively. The overlap is symmetric in its 
arguments h A and h B on ly if the signal and template, together with their parameters, are interchanged. That is, 
0[h A {\ A ), h B (X B )] = 0[h B {X B ), h A {\ A )] but, in general, 0[h A (X A ), h B {X B )} ^ 0[h A {X B ), h B (X A )}. Therefore, the 
maximized overlap e AB need not be symmetric. The process of maximization, in which the parameters of the "signal" 
are kept fixed and those of the "template" are varied, breaks down the symmetry. The lack of symmetry arises 
primarily because the signal manifolds M A . B representing the two families are distinct; the nearest "distance" from 
a coordinate point P on M. A to a point on M. B need not be the same as the nearest distance from P on M. B to a 
point on M. A . 

B. Effectulness, faithfulness and loss in event rates 

A direct measure of the efficiency of a template bank is the loss of event rates due to differences between the 
template family and the exact signal. The loss of event rates is determined by two factors: the effcctualness of the 
templates in matching the exact waveforms and the minimal match of the template bank itself. In this section, we 
will quantify this relation. 

In Fig.[2] 8 we sketch a portion of the waveform space. The solid line represents the template family subspacc. Dots 
represent various waveforms: (i) /iTmpitC^i) an d ^Tmpitfe) are two neighboring templates in the template bank with 
physical parameters x\ and X2] (ii) /iTmpit(so) and /iTmpit^o) are waveforms in the same family as the templates to 
be chosen as discussed below; (iii) h e (xo), /inr(^o) an d /ieob(^o) are exact, numerical and EOB waveforms of the 
same physical parameters xq, respectively. [The EOB waveform is calibrated to the numerical simulation.] We choose 
x' such that the overlap between /iTmpit (^l ) a nd /iTmpit(^o) is the minimal match (see below) of the template bank. 
We choose xq such that h e (xo) is the exact waveform that has larger overlap with /iTmpit (^o) than with any other 
waveforms in the template family. This overlap is larger than the one between h e (xo) and /iTmpit(^o) even though 
they have the same physical parameters, because of the systematic difference between the family of exact waveforms 
and the family of templates. 

We define the distance in the waveform space between two waveforms h and q by the scalar product -Jl — 0[h, q]. 
For convenience, we define the mismatch to be the square of the distance. The overlap between /iTmpit(^i) and 
^■Tmpit(^o) i s the minimal match and we denote the corresponding mismatch by dyiM = 1 — MM. Similarly, 1 — c?e 
and 1 — dp arc the effcctualness and faithfulness of the template family with the exact waveform h e (xo), respectively. 
The mismatch between h e (xo) and the closest template /iTmpit(^i) quantifies the reduction in signal-to- noise ratio 
when the template bank is used to search for the exact waveform. We denote this mismatch by dER. When these 
mismatches are small, by Pythagorean theorem, we have an the approximate relation djjR — c^mm + cfe- Assuming 
uniform spatial distribution of sources, the reduction in event rate is 1 — (1 — c?er) 3 — 3<iER. Therefore, if we want to 
satisfy the usual requirement of < 10% loss in event rate, we need g?er = ^mm + ^e < 3.5%. Typical minimal match 
adopted in current searches has either c^mm = 3% or g?mm = 1%, which means, in the first case, an extremely rigorous 
requirement on the effectualness: dE < 0.5%, or in the second case, a reasonable requirement of dE < 2.5%. The 
latter is achievable by PN models. Note that, if both the minimal match of a template bank and the effcctualness of 
the template model are 97%, the loss in event rate rises to 17%. 

However, it is not possible to calculate c?er since we do not know the exact waveform h e (xo). In this paper, we 
adopt two strategies to estimate g?er : (i) we calculate the mutual effectualness of PN models for low-mass binaries and 
assume it to be a good representation of their effectualness with exact waveforms; (ii) we approximate h e (xo) with the 
EOB waveform /ieob(^o) calibrated to the numerical simulations. We can verify the goodness of the latter assumption 
as follows. The mismatch between the best EOB waveforms [HI, |45[ and the numerical waveforms is less than 10~ 3 . 
In Ref. [45| , the authors calculated the mismatch among accurate numerical waveforms generated by simulations with 
different resolutions and/or extraction schemes, as well. They found that the mismatch is less than 10 -4 . We consider 
the latter as an estimate of the mismatch between exact and numerical waveforms. In the worst case, the mismatch 
between the exact and EOB waveforms with the same physical parameters is roug hly (Vl0- 3 + Vl0- 4 ) 2 = 1.7 x 10" 3 . 
Therefore, we can conclude that by approximating h e (xo) with /ieob(^o)- we underestimate the loss of event rate by 
at most 0.5%. 

Notice that the effectualness result presented in the following sections is slightly different from 1 — dg- It is 
obtained through discrete searches over template parameters using template banks with MM = 0.99 rather than 
through continuous searches. Therefore, the mismatch associated with this effectualness result includes already the 
discreteness effect in the template banks, i.e. a mismatch d^ M = 0.01. In this case, if a search is carried out with 



This figure is very similar to Fig. 3 of Ref. l2- r | 
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a template bank of a different minimal match, say MM= 1 — g?mm = 0.97, to calculate the loss of event-rate, a 
mismatch of g?mm — <^mm = 0.02, instead of 4m; needs to be added to the effectualness result in this paper, i.e. 
^er = d-MM — <^mm + ^E- The only exception in this paper is the effectualness result between EOB models presented 
in the Conclusions which is obtained through a continuous search. 



We have chosen three conventional systems, binary neutron stars (BNS), binary black holes (BBH) and binary neu- 
tron star-black hole systems, but we have chosen the BNS and BBH systems to be slightly asymmetric, (1.38, 1.42)M 
and (9.5, 10.5) M Q but NS-BH is chosen to be the conventional (10, 1.4) M Q system. To this we have added another 
binary with component masses (4.8, 5.2) which lies on the border line between where most PN families arc similar to 
one another and where they begin to differ. 

We compute overlaps maximized over a template bank between seven different models (Taylor TT, TaylorT2, 
TaylorT3, TaylorT4, TaylorF2, TaylorEt, EOB), each at three different PN orders (v 4 , v 6 , v 7 ). The results will be 
presented in the form of a set of Figures. For each mass pair there will be one Figure consisting of 9 panels (one 
panel for each PN order), each panel containing seven curves (one each for each template family at that order) and 
each curve with 21 data points corresponding to signals from the seven PN families at each three different PN orders, 



We will present the results of our investigation in two complementary ways. We will first discuss the effectualness 
of the different PN families with each other. Such an analysis will help us understand how well the PN approximation 
has "converged" for the selection of detection templates. We then go on to look at the effectualness of the different 
approximants with the EOB signal that contains not only the inspiral but also the merger and ringdown parts. The goal 
of the latter analysis is to identify the region in the parameter space where one can safely use any PN approximant 
template in a search, without worrying about the loss in signal-to-noise ratio that might arise due to our lack of 
knowledge of the real signal, but without expending undue computational resources. Outside this region, however, 
one must use template families that are calibrated to waveforms obtained from numerical relativity simulations. 



The effectualness of the different PN families with each other is shown in Figs. [3] (Initial LIGO) and 2] (Advanced 
LIGO) for four different systems with component masses as indicated at the top of each sub-figure. In each sub-figure, 
the top panels correspond to the effectualness of different template families at 3.5PN order, middle panels to 3PN 
order and bottom panels to 2PN order. For each template family considered we find their overlap with signals from 
different PN orders (as indicated along the x-axis) and approximants (as indicated by the text Tl, T2, etc.). Each 
symbol corresponds to the overlap obtained by a different template family: (black) circles to TaylorTl, (red) squares 
to TaylorT2, etc., with signals from different PN families. Note that we have used the logit scale 9 for the vertical axis. 
This is so that (minor) disagreements between the different approximants are made clearly visible. Note that since 
we are considering systems with low total mass, say < 20M Q , in this section we use the EOB model terminated at 
the EOB light ring, that is we do not include the merger and ringdown parts. 

Conventionally, one says that two approximants A and B are in close agreement with each other if their mutual 
effectualness 6ab is 0.965 or greater [171 ] . Since in this study we are dealing with a rather large number of different 
PN families (21 in all), we shall relax this condition a bit to 0.95. However, we shall indicate in Sec. I VIII the 
region of the parameter space where the effectualness is better than 0.965, but we shall also quote regions where 
the effectualness drops to a low value of 0.9. The latter should be helpful for data analysis pipelines that employ a 
multi-stage hierarchical search, the first stage of which deploys a coarse grid of templates. 

These figures reveal many different aspects of the (dis) agreements between the different approximants but we shall 
only mention in our discussion the "diagonal" behaviour, i.e. overlaps of each template family with a signal family 



Recall logit (p) = log (jz^J • 




C. Choice of binary systems and PN orders 



2PN, 3PN and 3.5PN. 



VI. RESULTS OF THE EFFECTUALNESS OF PN TEMPLATES 



A. Mutual effectualness of various PN Inspiral template banks 
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FIG. 3: The plot shows the effectualness of templates and signals of different post-Newtonian familes and orders for four 
different binary systems for Initial LIGO. For a template from a given PN approximation (indicated by different line styles and 
symbols) and order (top panel 3.5PN, middle panel 3PN and bottom panel 2PN) we compute the effectualness of each of the 
templates with signals from each of the seven families, TaylorTl (Tl), TaylorT2 (T2), Taylor T3 (T3), Taylor T4 (T4), TaylorF2 
(F2), TaylorEt (Et) and Effective-One-Body (EOB), at 2PN, 3PN and 3.5PN orders. For instance, solid lines with filled circles 
give the effectualness of TaylorTl templates at 3.5PN (top panel), 3PN (middle panel) and 2PN(bottom panel) PN orders, 
with signals that belong to different PN approximations and orders. In clockwise order the panels from top left correspond to 
binaries consisting of two neutron stars, with masses 1.38 Mq and 1.42 Mq, two black holes with masses 4.8Mq and 5.2 Mq, 
two black holes with masses 9.5 Mq and 10.5 Mq and, finally, a neutron star and a black hole binary with component masses 
1.4 M and 10 Mq. 



from the same PN order. Focusing first on the Initial LIGO results (Fig. [3]), we see the evidence for the clustering 
of the various approximants at 3PN and 3.5PN orders for systems with a smaller total mass. In the case of BNS 
with component masses (1.38, 1.42) Mq, 2PN "diagonal" overlaps are dispersed between 0.74 to 1, 3PN and 3.5PN 
overlaps are all above 0.95, with TaylorEt having the smallest overlaps. 

In the case of BBH with component masses (4.8, 5.2) Mq, 2PN overlaps are between 0.8 and 1, 3PN overlaps are all 
greater than 0.95 except TaylorEt, 3.5PN overlaps are greater than 0.95 for all except TaylorEt, TaylorT3 and EOB. 
There are several important points to note: As discussed in Sec. IIV1 TaylorT3 terminates somewhat prematurely 
before reaching the last stable orbit. Therefore, one expects to have poorer overlaps for all templates if TaylorT3 
signal terminates in band, which will be the case for systems with a total mass greater than about 10 Mq. The 
asymmetry in the overlaps mentioned in Sec. IV Al is apparent in the case of TaylorEt: The overlaps of all templates 
with TaylorEt signal is greater than the converse, namely the overlaps of the TaylorEt templates with other signals. 
The poorer performance of EOB templates (terminated at the light ring) is due to the fact that the waveform has 
power in band beyond the last stable orbit. 

In the case of NSBH with component masses (1.4, 10) Mq, 2PN "diagonal" overlaps are distributed between 0.6 
and 1, 3PN and 3.5PN overlaps are consistently above 0.95 except for TaylorEt signals (both orders) and TaylorT3 
(at 3.5PN). 

In the case of BBH with component masses (9.5, 10.5) Mq, there is no agreement between approximants irrespective 
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FIG. 4: Same as Fig. but for Advanced LIGO. 



of the PN order. In this sense, one cannot trust using any particular approximant as a. search template. 

Let us now turn to Fig. 2] which depicts the results for Advanced LIGO noise power spectral density. In the case 
of BNS with component masses (1.38, 1.42) Mq, the 2PN "diagonal" overlaps are between 0.4 and 1 (note that some 
of the data points are below the scale of 0.5 that we employ). The 3PN (except TaylorEt signal) and 3.5PN (except 
TaylorT3 template and TaylorT3 and TaylorEt signals) overlaps are uniformly larger than 0.95. The effcctualncss 
of all templates with TaylorEt signal is generally smaller (0.6-0.8) than the effectualness with a TaylorEt template. 
In the case of BBH with component masses (4.8, 5.2) Mq, the 2PN overlaps could be as small as 0.65. At 3PN, all 
approximants (except TaylorEt templates) and 3.5PN (except TaylorEt and TaylorT3 templates) the overlaps are 
0.95 or greater. In the case of NSBH with component masses (1.4, 10) Mq, the 2PN overlaps are as low as 0.4. At 
3PN and 3.5PN, the overlaps are larger than 0.95 except in the case of TaylorEt signals (3PN, 3.5PN) and TaylorT3 
templates (3.5PN). In the case of BBH with component masses (9.5, 10.5) Mq, the 2PN overlap could be as low as 
0.7. The overlaps are larger than 0.95 at 3PN except in the case of EOB templates and TaylorEt and EOB signals. 
Finally, at 3.5PN order the different approximants are seen not to agree with each other very well. The cause of these 
features is the same as our discussion for Initial LIGO. 



B. Discussion 



In the case of binary neutron stars, the merger occurs far outside the sensitive band of the detector and even the 
late stages of inspiral is out of band. Binary neutron stars will very much be in the adiabatic regime as the signal 
sweeps through the band and a good test of the PN approximation is to ask how well the different waveforms agree 
with one another in this regime. The finite bandwidth of the detector essentially probes this regime for binary neutron 
stars. Note that the effcctualncss amongst different PN families at 2PN order is pretty poor but greater than 0.95 
(with the exceptions discussed earlier) at 3PN and 3.5PN orders. In the case of Advanced LIGO (cf. Fig. QJ, the lower 
frequency cutoff used in computing the overlap integrals is 20 Hz and a binary neutron star spends more than 750 



18 



I 
o 



1 41.38, 1.42) M 



O/*-* * * * * 

T1T2 T3T4 F2 Et EoB 










1 






-10 




< 






o 






c 


-20 








-(4.8, 5.2) M Q 




-30 



■ 
-10 ■ 



O -20 - 




-30 L 



-(9.5, 10.5) M 



20 


-20 h (1.4, 10) M, 



1 - 

l 



-1 

30 

> 20 

O 10 
o 


30 



> 20 

> 

< 

O 10 
o 





20 

lo 



-(1.38, 1.42) M 



Tl 
T2 
T3 
T4 
F2 
T-TEt 
►-►EoB 





• 


X4.8, 5.2) M 






A 

/ ^ 








(TT- 



(9.5, 10.5) M, 




_(1.4, 10) M Q 





0.5 




10 









-10 


< 


c 




c 


-20 






-30 




10 









-10 


< 


c 




c 


-20 






-30 




20 




10 







o 




c 


-10 






-20 



(1.38, 1.42) M Q 
Tl T2 T3 T4 F2 Et EoB 



. (4.8, 5.2) M Q 




7(4.8, 5.2) M Q 





(9.5, 10.5) M n 



(1.4, 10) M Q 



0.5 

> 
5 

< o 

o 

o 



-0.5 




30 




20 


> 






10 


o 




o 









10 




30 


> 


20 


1 






10 


o 




o 









10 




20 


> 





1 






20 


o 




o 


40 






60 



(1.38, 1.42) M Q 



(4.8, 5.2) M, 




v. 



(9.5, 10.5) M n 



w -w > 

f (1.4, 10) M. 



Tl 
■< <■ T2 

♦■-♦T3 

<-< F2 
T-TEt 
►-►EoB 



FIG. 5: Percentage bias in the estimation of the total mass M and symmetric mass ratio v at 3.5PN order. Left two columns 
are for Initial LIGO and the right two for Advanced LIGO. The bias AM/M is defined as AM = (1 - Mrmpit/M S g„i) , where 
A/sgni and MTmpit denote the total mass corresponding to the signal and the template that obtained the maximum effectualness, 
respectively (and similarly for v). What is plotted is percentage bias. The bias arises because the template family (as indicated 
in the key) is different from that of the signal family (as indicated in the top left panel as Tl, T2, etc.). 



cycles in band. Effectualness of 0.95 or greater means that the waveforms remain in phase over the entire duration of 
the signal. Of course, in reality the parameters of the signal and the template are not the same, but even so this is a 
remarkable success of the PN scheme. 

For a BBH system with masses (4.8, 5.2) M©, we see that 2PN and 3PN order templates are qualitatively similar 
to the binary neutron star case. However, we can see a marked deterioration of the effectualness at 3.5PN order. For 
a system of total mass of 10 M©, the Schwarzschild LSO occurs at ~ 440 Hz and the detector is sensitive to the late 
stages of the inspiral phase. It is not entirely surprising, therefore, that different PN orders do not agree with each 
other to the same extent as in the binary neutron star case. However, note that, with the exception of TaylorT3, 
which terminates at a frequency somewhat lower than others, and TaylorEt, all other templates have effectualness of 
0.95 or better with each other. Among approximants that agree with each other, EOB has the smallest effectualness. 
This is because the latter model contains the plunge phase of the coalescence with ending frequencies far higher than 
the LSO while other approximants do not have the plunge phase. 

The LSO of a BBH with component masses (9.5, 10.5) M©, is at ~ 220 Hz and the plunge phase spans 220 Hz to 
about 600 Hz. Therefore, the detector is pretty sensitive to the late phases of the coalescence. We see deterioration of 
the effectualness, both at 3PN and 3.5PN orders. Apart from TaylorT3, whose poor overlaps at 3.5PN are explained 
by the early termination of the signal, the EOB stands out by achieving overlaps as low as 0.92 with other families. 

As a final example, the effectualness of templates for a signal from a neutron star-black hole binary of masses 
(1.4, 10) Mq, we see that the different PN families, including the EOB, are in good agreement with each other, with 
the sole exception of TaylorEt. In fact, the convergence amongst different families seems to be somewhat better than 
the BBH system of component masses (9.5, 10.5) Mq. 

At this juncture, it is worth pointing out that our numerical results for effectualness in the subset of cases where 
TaylorEt is chosen as the signal model, are consistent with those in Ref. [4Vj j . which investigated the fitting factors to 
ascertain if 3.5PN TaylorEt signals could be effectually and faithfully searched by TaylorTl, TaylorT4 and TaylorF2 
templates. There is agreement too on the general features of our results with regard to systematic biases, the 
dependence on the total mass and qualitative factors underlying them. However, this agreement of numerical results 
for faithfulness and effectualness in no way extends to the general motivation and claims regarding the TaylorEt 
approximants (47l . |49| . [50j and, hence, are worth clarifying. 

Indeed, there is no basis to refer to the x-bascd orbital phasing equation Eq. (|3.5ap as Newtonian [5(|, since 
the u> here is nPN accurate (depending on the PN-generation order one is working at) and implicitly incorporates 
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FIG. 6: Overlaps of different 3.5PN approximants with the EOB inspiral-merger-ringdown signal in Initial LIGO in the 
(mi — m,2)MQ plane. The approximants considered from left-to-right are TaylorTl, TaylorT2, Taylor T3 (top panels), and 
Taylor T4, TaylorF2, TaylorEt (bottom panels). The contours correspond to overlaps of 0.9, 0.95 and 0.965. 



conservative contributions to gravitational-wave phase evolution at various PN orders. It is incorrect to claim [47] 
that conservative contributions to the gravitational-wave phase evolution do not appear in the standard approximants, 
or that the TaylorEt-bascd scheme treats conservative and radiation-reaction contributions more equitably than the 
standard x-based approximants. It is misleading [47l | to refer to only TaylorEt-based approximants as "fully gauge 
invariant in contrast to EOB" (especially in the circular orbit case). All cc-based schemes are also fully gauge invariant. 
Finally, one may work in specific convenient coordinate systems as do EOB and numerical relativity simulations, as 
long as one deals with and compares gauge invariant quantities at the end. 

In our view, the very different behaviour of the TaylorEt approximant relative to the standard x-based approximants 
may be traced to the manner in which the orbital phasing is "packaged" in the two schemes. In the x-based schemes the 
orbital phasing is implicitly in a resummed form, since the phasing is written in an appropriate PN-accuratc angular 
velocity w„pn (n = 2, 3 for 2PN, 3PN templates). On the other hand, the representation in terms of £, relative to the 
x schemes, is a re-expanded form. And indeed, based on the comparison between analytical schemes and numerical 
relativity simulations, the £ schemes do relatively worse. The feature related to the monotonic-convergence of the 
TaylorEt scheme is of secondary importance in comparison to the main requirement of high phasing accuracy of an 
analytical model with numerical relativity simulations over all mass-ratios. 

A few general comments are in order before we conclude this Section. We do not believe that at present there are 
convincing theoretical reasons to consider any one particular PN family of inspiral models to be a privileged signal 
model. Consequently, the best that one can do is to examine the mutual closeness of these various inspiral models, 
as we have done, and work at the PN order where these various template families display the greatest agreement. It 
is precisely in this regard that the viewpoint we present here differs from those in (47l. I49t[5f3| which assumes primacy 
for one specific approximant, namely the TaylorEt approximant, based on theoretical motivations that at present do 
not appear to be fundamental or compelling. Consequently, though there is no difference in the numerical results in 
the subset of cases that are common in our investigations, there is a big difference in the conclusions that we believe 
can be inferred. For instance, before one can legitimately decide on the inability of standard template banks in the 
gravitational data pipeline to detect signals from binaries with eccentricity [50], it is necessary to first fold in the 
differences in the simpler quasi-circular case arising on account of different parametrisations. Similar considerations 
should be borne in mind when dealing with analogous problems in the spinning case. 

Based on the analysis presented heretofore, we conclude that the PN approximation has pretty much converged at 
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3PN and 3.5PN orders 10 , as long as the total mass is less than about 12 M Q (with the exceptions discussed in the 
previous Section). 

For heavier binaries, the approximants begin to differ considerably, and this is almost entirely because the adiabatic 
approximation begins to breakdown and the plunge and the merger phases become more and more important. Hence, 
in the next Section we will supplement the present analysis by looking more precisely into the overlaps of the different 
PN templates with a prototype of the more complete signal model, namely the EOB model, including the merger and 
ringdown parts. 



C. Biases in the estimation of parameters 



Recall that, in the computation of the effcctualness one maximizes the scalar product of a (normalized) signal with 
a template over the parameters of the template keeping those of the signal fixed. Therefore, one can get an idea of 
how dissimilar the parameters of an approximant need to be in order to match a given signal. This is a systematic 
effect that leads to a bias in the estimation of parameters if the template approximant is not the same as the signal 
approximant. Let the total mass of the signal and template waveforms be, respectively, Ms gn i and Mxmpit, when the 
scalar product is maximized. The percentage bias AM in the total mass is defined as AM = 100(1 — Mxmpit/Msgni), 
and similarly for the symmetric mass ratio v. 

For a given binary, the biases are qualitatively similar for Initial and Advanced LIGO noise power spectral densities. 
In general, the biases are appreciably smaller at 3PN and 3.5PN orders than at 2PN order and progressively increase 
with the total mass, although they are far larger than the statistical errors computed using the Fisher information 
matrix [23|. Figs. [5] plots the percentage biases in the total mass M and symmetric mass ratio v at 3.5PN order. 
The left two (right two) columns use the Intial LIGO (Advanced LIGO) noise spectral density. For the four systems 
considered, namely (1.38 , 1.42) M Q , (4.5 , 5.2) M , (1.4 , 10) M , and (9.5 , 10.5) M Q binaries, the largest bias in the 
total mass M is 1%, 20%, 20% and 20%, respectively, and the symmetric mass ratio v is 1%, 25%, 70% and 25%, 
respectively. 



Though qualitatively we may expect similar results for Virgo, quantification requires an analysis using the Virgo noise curves. Needless 
to add, that the situation for a space detector like LISA can be expected to be even more different and interesting to study 
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VII. RESULTS OF THE EFFECTUALNESS OF PN TEMPLATES WITH THE FULL WAVEFORM 

Having established the convergence of PN approximations at 3PN and 3.5PN orders (for determining effectual 
templates for detection) in the regime where the approximation is expected to be valid, let us now examine the region 
in the parameter space where PN families can be used as search templates. To achieve this goal we will use the EOB 
model calibrated to numerical relativity simulations [39| . For brevity, we have omitted plots of the effectualness of 
the 3PN approximants with this EOB model; they are quite similar to the 3.5PN plots. 

Although Ref. [3t| explored the agreement between the EOB model and numerical simulations for several modes, 
in this study we will work with only the dominant harmonic (i.e., the mode) at leading PN order. Higher-order 



amplitude corrections are known to be important for parameter estimation [60l . l6l| and a future study must repeat 
this investigation with the full waveforms. 

Fig. [5] shows the effectualness of the six PN families TaylorTl, TaylorT2, TaylorT3 (top panels, respectively from 
left to right), TaylorT4, TaylorF2, and TaylorEt (bottom panels, respectively from left to right) for Initial LIGO noise 
power spectral density. Fig. [7] shows the same but for Advanced LIGO noise power spectral density. The effectualness 
was computed using a hexagonal template bank [7^ | and is shown as a gray-scale map in the space of the component 
masses that are taken to vary from 3Af Q to 14.5 M Q . For all the maps we have chosen the gray-scale to vary from 
0.76 to 1. The dotted contours show effectualness at three values: 0.965, 0.95 and 0.90. 

The trends of the overlaps is rather similar irrespective of which noise power spectral density we use, although the 
actual overlaps are systematically smaller in the case of Advanced LIGO as compared to Initial LIGO. This is due 
to the broader frequency sensitivity of the former in relation to the latter. The following discussion is, therefore, 
applicable in both cases. 

Let us first note some peculiarities. TaylorT3 at 3.5PN leads to particularly ineffectual templates. As mentioned 
before, TaylorT3 at 3.5PN terminates rather prematurely. The LSO defined by the Schwarzschild potential is at 
/lso ~ (440/10 Mq) Hz, but TaylorT3 at 3.5PN approximants terminate at ~ (220/10 M Q ) Hz. This discrepancy is 
so large that even with the biases in the component masses allowed in the computation of the effectualness (recall that 
we maximize the overlap over template masses), which, in principle, makes it possible for a template of a lower mass 
to match a signal of a higher mass, TaylorT3 is unable to achieve good overlaps. This is because a mismatch in the 
component masses can make a template more, or less, asymmetric than the signal, which has the effect of increasing, 
or decreasing, the duration of the template relative to the signal. While small differences in the ending frequencies 
can be achieved by a mismatch in the total mass without affecting the signal duration too greatly, large differences 
cannot be compensated by such a mismatch in the parameters. 

At 3PN and 3.5PN the effectualness of TaylorEt with a EOB signal for a binary of component masses (3, 10) M Q 
[respectively, (10, 10) M Q ] is 0.83 and 0.90 [respectively, 0.87 and 0.89]. This is because amongst all PN approximants 
TaylorEt seems to converge far slower than any other. Further, an examination of the coefficients in the PN terms 
of the phasing formulas in Eqs. (|3.13p indicates that higher order PN terms have increasingly greater coefficients. In 
general, it has been observed that the appearance of such larger coefficients in higher order terms of an approximant 
scheme inevitably worsens its convergence and the present instance may be no exception to this case 11 . 

With the exception of the peculiarities noted above, we see that all approximants do progressively better at higher 
PN orders. Conclusions drawn in the previous Section with regard to the convergence of the PN approximations are 
further corroborated here where we have measured the overlaps with a signal that is matched to numerical relativity 
simulation, which can, therefore, be taken to be close to what a real signal might be. 

Computationally, TaylorF2, with its phasing formulas given explicitly in the Fourier domain, is the least expensive. 
This is because matched filtering is most easily carried out in the Fourier domain, which means that a time-domain 
approximant must be Fourier transformed before computing the cross correlation. By employing TaylorF2 models 
one can avoid one forward Fourier transform. Moreover, TaylorF2 offers the flexibility in the choice of the ending 
frequency. Unlike the time-domain models, which have either a natural ending frequency defined by the extremum of 
the binding energy or the frequency evolution stops before reaching LSO, TaylorF2 has no such restriction. In fact, 
as obtained in Refs. [H, Hlj|, by extending the upper cutoff beyond the usual upper cutoff (i.e., the Schwarzschild 
LSO), the TaylorF2 model matches remarkably well with numerical relativity waveforms for a far greater range of 
masses. However, as noted in Rcf. 36] the ending frequency that must be employed in order to achieve the best 
match with numerical-relativity waveforms depends on the noise power spectral density. This could turn out to be an 
unnecessary computational burden in a data analysis pipeline. The alternative is to choose the upper frequency cutoff 
as an additional search parameter or allow unphysical values of v > 0.25 [H, [H, |47j or to include a p4PN term in the 



11 While comparing the coefficients it may be useful to note that v ~ l/y/E corresponds to f in the range of 0.136-0.138 depending on the 
symmetric mass ratio v and the PN order. 
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template phase and calibrate it to numerical simulations (35| . The first two choices would result in an unwarranted 
increase in the computational cost of a search as also in the false alarm rate, and we advice against it. The third 
choice could be pursued, but it should be augmented by a more complete description of the merger/ringdown signal 
— for example by introducing a slope break in the waveform amplitude and a superposition of Lorentzians [34l . [35| 

If a search requires the minimal match to be much smaller than 0.95 (as, for example, in a hierarchical search) one 
can extend a search with TaylorF2 to a total mass of 20 Mq with effectualness of 0.90. 

Before advanced detectors begin to operate, there will be a period when LIGO and Virgo will operate with sensitivi- 
ties slightly larger than, but bandwidths similar to initial detectors - the so-called Enhanced LIGO and Virgo+. Since 
Virgo and Virgo+ are expected to have a sensitivity bandwidth similar to Advanced LIGO the results presented in 
this paper are qualitatively similar to in those cases too. Moreover, as our results are only sensitive to the bandwidth, 
conclusions drawn by using the noise spectral density of Initial LIGO will also be valid for Enhanced LIGO. 

All approximants (no exceptions) achieve an effectualness of 0.95 or better at 3PN and 3.5PN orders, for binaries 
whose total mass is less than about ~ 12 Mq. From the view point of effectualness alone, we conclude that searches 
for binary black holes, in Initial, Enhanced and Advanced LIGO, could employ any of the 3PN or 3.5PN families as 
long as the total mass is smaller than about ~ 12 Mq. The final choice of the PN family should be based on other 
criteria. If it is desired that the minimal match of a template bank is 0.965 or greater, then the best strategy would 
be to use the full EOB waveform calibrated to numerical relativity. 

Another criteria to be considered is the computational cost. A typical matched filter search in LIGO data must 
compute thousands of template signals for every 2048 second data segment. This can be a heavy burden if it takes 
a significant amount of time to compute each template. The EOB templates are computed in the time domain by 
solving a set of differential equations, and the frequency domain signal is then computed via Fourier transform. For 
low-mass systems this cost can become significant and will of course vary depending on the implementation and 
hardware used. 

We have estimated the cost to compute TaylorF2 and EOB templates using their implementation in the LIGO 
Algorithm Libraries (LAL) code used for matched filtering searches in LIGO data. We find that for a total mass 

> 40 Mq, the EOB templates take a factor of 2 longer to generate than the same TaylorF2 signals. For a (10, 10) Mq, 
(5, 5) Mq and (1.4, 1.4) Mq binary, the EOB templates take about a factor of 3, 7 and 20, longer to generate, 
respectively. We tested the waveform generation on a high performance computer with 32 2.7 GHz CPUs and 132 
GB of RAM. On this system, EOB templates with a total mass > 40 Mq can be generated in about 0.1s, while the 
(10, 10) Mq EOB template could be generated in about 0.5s. Since LIGO searches employ thousands of CPUs, this is 
feasible. However, for lower mass signals, the time needed grows rather quickly and about 4s are needed to compute the 
(1.4, 1.4) Mq EOB template. It may be possible to reduce the computational cost somewhat by optimizing the EOB 
waveform generation code, but the lowest mass templates would almost certainly still have a significant computational 
cost. Thus, the increased computational cost must be weighed against the benefit of increased effectualness for lower 
mass signals. 

VIII. FAITHFULNESS 

For completeness, we also report on the faithfulness of the different PN approximants with respect to one another. 
The faithfulness is the overlap between normalized template and signal approximants when maximizing only over 
the time and phase at coalescence, tc and 4>c- In Tables HT1 and ITTTl we list the faithfulness for each pair of PN 
approximants at their highest PN order, that is 3.5PN order, except for the EOB model which uses a p4PN order 
coefficient, for both Initial and Advanced LIGO and for each of our reference binaries. 

In the first row and column of the left panel of Table QT] notice that every approximant has an overlap of at least 
0.97 with the EOB model for both Initial and Advanced LIGO. That all approximants have good agreement for a low 
mass binary without searching over mass parameters is further evidence that the 3.5PN approximants are rather close 
to one another during the adiabatic inspiral. Note that the T2, T3, T4 and F2 approximants all have a faithfulness 

> 0.99 with the EOB model, while the Tl and Et approximants have somewhat worse agreement at about 0.97. For 
each pair, the faithfulness for Initial and Advanced LIGO are quite similar for these low mass binaries. 

In the right panel of Table [Til we increase the total mass to 10M Q while keeping the mass ratio nearly equal. The 
faithfulness drops for every pair of approximants as the merger begins to enter the sensitive band. Recall that for these 
masses, all pairs of approximants can achieve an effectualness of at least 0.95 by searching over the mass parameters. 
When we fix the masses, the T2, T4 and F2 approximants still have very good agreement with the EOB model, 
with faithfulness of 0.97-0.98. The EOB-T3 faithfulness has degraded somewhat to 0.93 - 0.94, and the Et and Tl 
approximants have rather poor agreement with the EOB model with faithfulness in the range 0.84 — 0.92. Note that 
the faithfulness is typically lower for Advanced LIGO than for Initial LIGO. We attribute this to the signals having 
a longer duration (and thus more time to accumulate a phase difference) in Advanced LIGO's wider sensitivity band. 
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TABLE II: Faithfulness of different approximants for (1.42, 1.38) M Q (left panel) and (5.2, 4.8) M Q (right panel) binaries. The 
rows label template approximant, while the columns label signal approximant. For each pair, the top number is Initial LIGO 
while the bottom number is Advanced LIGO. All approximants are at 3.5PN order, except our EOB model which has a p4PN 
coefficient. 





EOB 


Tl 


T2 


T3 


T4 


Et 


F2 




EOB 


Tl 


T2 


T3 


T4 


Et 


F2 


EOB 


1 


.969 


.994 


.997 


.990 


.970 


.994 


EOB 


1 


.916 


.974 


.938 


.981 


.888 


.970 




1 


.971 


.996 


.998 


.991 


.974 


.996 




1 


.877 


.973 


.928 


.978 


.841 


.968 


Tl 


.969 


1 


.982 


.981 


.987 


.928 


.982 


Tl 


.916 


1 


.974 


.926 


.964 


.784 


.975 




.971 


1 


.984 


.983 


.990 


.920 


.984 




.877 


1 


.955 


.892 


.947 


.653 


.957 


T2 


.994 


.982 


1 


.998 


.999 


.958 


1.000 


T2 


.974 


.974 


1 


.949 


.993 


.861 


.993 




.996 


.984 


1 


.999 


.999 


.961 


1.000 




.973 


.955 


1 


.932 


.994 


.775 


.995 


T3 


.997 


.981 


.998 


1 


.997 


.959 


.998 


T3 


.938 


.926 


.949 


1 


.943 


.925 


.944 




.998 


.983 


.999 


1 


.998 


.961 


.999 




.928 


.892 


.932 


1 


.926 


.876 


.930 


T4 


.990 


.987 


.999 


.997 


1 


.950 


.999 


T4 


.981 


.963 


.993 


.943 


1 


.854 


.995 




.991 


.990 


.999 


.998 


1 


.949 


.999 




.978 


.947 


.994 


.926 


1 


.766 


.996 


Et 


.970 


.928 


.958 


.959 


.950 


1 


.958 


Et 


.888 


.785 


.861 


.925 


.854 


1 


.852 




.974 


.920 


.961 


.961 


.949 


1 


.961 




.841 


.653 


.775 


.876 


.767 


1 


.770 


F2 


.994 


.982 


1.000 


.998 


.999 


.958 


1 


F2 


.970 


.975 


.993 


.944 


.995 


.853 


1 




.996 


.984 


1.000 


.999 


.999 


.961 


1 




.968 


.957 


.995 


.930 


.996 


.770 


1 



TABLE III: Same as Table UU but for (10.5, 9.5) M© (left panel) and (10, 1.4) M (right panel) binaries. 





EOB 


Tl 


T2 


T3 


T4 


Et 


F2 




EOB 


Tl 


T2 


T3 


T4 


Et 


F2 


EOB 


1 


.877 


.882 


.650 


.923 


.860 


.910 


EOB 


1 


.977 


.973 


.817 


.859 


.526 


.990 




1 


.811 


.864 


.721 


.910 


.775 


.889 




1 


.959 


.972 


.801 


.797 


.413 


.993 


Tl 


.877 


1 


.972 


.712 


.970 


.817 


.982 


Tl 


.977 


1 


.972 


.796 


.805 


.508 


.991 




.811 


1 


.955 


.785 


.943 


.638 


.966 




.959 


1 


.954 


.753 


.691 


.398 


.978 


T2 


.882 


.972 


1 


.742 


.968 


.886 


.959 


T2 


.973 


.972 


1 


.835 


.894 


.543 


.980 




.864 


.955 


1 


.831 


.969 


.784 


.959 




.972 


.954 


1 


.820 


.834 


.430 


.976 


T3 


.650 


.712 


.742 


1 


.707 


.716 


.709 


T3 


.817 


.796 


.835 


1 


.851 


.778 


.818 




.721 


.785 


.831 


1 


.794 


.782 


.790 




.801 


.753 


.820 


1 


.841 


.631 


.798 


T4 


.923 


.971 


.968 


.707 


1 


.906 


.986 


T4 


.859 


.805 


.894 


.851 


1 


.595 


.852 




.910 


.943 


.970 


.794 


1 


.785 


.988 




.797 


.691 


.834 


.841 


1 


.456 


.779 


Et 


.859 


.817 


.886 


.716 


.906 


1 


.845 


Et 


.526 


.508 


.543 


.778 


.595 


1 


.525 




.776 


.639 


.784 


.783 


.785 


1 


.707 




.413 


.398 


.430 


.631 


.456 


1 


.411 


F2 


.909 


.982 


.959 


.708 


.985 


.846 


1 


F2 


.990 


.991 


.980 


.818 


.852 


.525 


1 




.889 


.967 


.959 


.790 


.988 


.706 


1 




.993 


.978 


.976 


.799 


.779 


.411 


1 



In the left panel of Tabic ITTTl we increase the total mass to 2QMq while again keeping the mass ratio nearly equal. 
Once again, the faithfulness drops for all cases as the merger and ringdown become more important. The T4 and F2 
approximants have the best agreement with EOB, they are the only approximants to achieve an overlap greater than 
0.9 with EOB. The overlap between T3 and EOB has dropped dramatically to 0.65 and 0.72 for Initial and Advanced 
LIGO respectively. 

The right panel of Tabic ITTTl gives the faithfulness for each approximant pair for an asymmetric (10, 1.4) Mq binary. 
The EOB-F2 faithfulness is very good at 0.99. The Tl and T2 approximants also have good agreement with the 
EOB model with faithfulness 0.96 — 0.98. The T3 and T4 have poor agreement with the EOB model with faithfulness 
0.80 — 0.86. For this mass pair, the Et approximant has very poor agreement with all of the others, the faithfulness 
is < 0.60 for every approximant except T3. 

We see a clear trend of decreasing faithfulness as the total mass of the binary increases. This is due to the late 
inspiral, merger and ringdown moving into the sensitive band and becoming more important for higher mass binaries. 
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FIG. 8: Effectualness (left y-axis) and the corresponding loss in event rate (right y-axis) of 3.5PN approximants with the EOB 
inspiral-merger-ringdown signal calibrated to numerical relativity in Initial LIGO (left panel) and Advanced LIGO (right panel) 
as a function of total mass for 1:1, 4:1 and 10:1.4 mass ratios. The EOB curve is the effectualness between the uncalibrated 
3.5PN EOB model containing only the inspiral and the calibrated inspiral-merger-ringdown EOB signal. Note that any of these 
approximants are suitable for detection templates below a total mass of about 12 Mq for both Initial LIGO and Advanced 
LIGO, provided a 10% loss of event rate is deemed acceptable. 

The faithfulness is typically lower for Advanced LIGO than Initial LIGO due to its broader sensitive band. The 
faithfulness can vary with mass ratio. For example, for the (10, 1.4) Mq binary, the Tl and T2 approximants have a 
better faithfulness with the EOB model than the T4 approximant. However, for the nearly equal mass binaries, the 
T4 approximant has the greater faithfulness with the EOB model. The TaylorF2 approximant is generally the most 
faithful to the EOB aproximant, with one of the highest overlaps in each case. This is another argument for using 
TaylorF2 templates in the mass regime where EOB templates are too computationally expensive to be employed in 
a matched filtering search. 

IX. CONCLUSIONS 

In this paper we have examined the convergence of the PN approximation with the view to validating their use in the 
search for compact binaries in Initial, Enhanced and Advanced LIGO. We considered seven different approximants, 
each at three different PN orders, a total of 21 waveforms in all. We computed the effectualness of each of the 
waveforms with every other at 2PN, 3PN and 3.5PN orders by using a template bank constructed with a minimal 
match of 0.99 and Initial and Advanced LIGO noise power spectral densities. Our results from a sample of four 
binaries show that different PN approximations are consistent with one another at 3PN and 3.5PN order. They begin 
to differ only when the mass becomes so large that the plunge phase, not contained in standard PN waveforms in the 
adiabatic approximation, enters the detector band. 

The above conclusion is best summarized by Fig. [8l where we plot the effectualness of the various PN approximants 
(except for TaylorT3 and TaylorEt that we recommend be discarded, since we have shown that not only do they differ 
considerably from the others but importantly have poorer overlaps with EOBNR waveforms) with an EOB inspiral- 
merger-ringdown signal as a function of the total mass of the binary. These plots are convenient for identifying the 
-^crit above which the PN approximants begin to differ with one another. We find that any of the above approximants 
could be used as detection templates with less than a 10% loss in event rate up to a total mass of 12A/ Q for both 
Initial and Advanced LIGO. Note that this value of M cr i t is limited by the equal-mass case, as the value of M cr it 
corresponding to a 10% loss in event rate is somewhat greater for mass ratios of 4:1 and 10:1.4. We attribute this 
result to asymmetric binaries accumulating more signal at low frequencies than in the equal-mass case. Thus, for a 
fixed total mass, the merger and ringdown are less important for asymmetric binaries than for equal-mass binaries. 
Therefore, we conclude that we can safely use any of the above 3.5PN families as search templates to detect binaries 
whose total mass is less than about 12 Mq. However, purely from the point of view of computational burden TaylorF2 
is the least expensive and we recommend that TaylorF2 at 3.5 PN order be deployed as search templates below a 
total mass of 12 Mq. It is quite remarkable to note that up to a total mass of 30M Q , the uncalibrated EOB model at 
3.5PN order is rather close to the calibrated EOB inspiral-merger-ringdown signal. In fact, Ref. [43| found a phase 
difference of only 0.05 rads after 30 GW cycles, at roughly 3 GW cycles before merger between the EOB at 3.5PN 
order and the highly accurate equal-mass numerical waveform of Caltech/Cornell collaboration. 

For systems with total mass larger than about 12 Mq, TaylorF2 at 3.5PN might be effectual if the upper cutoff 
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TABLE IV: Effectualness and faithfulness of the EOB fiducial model [3S|| used in this paper (and currently employed by Initial 
LIGO) to the most recently improved EOB model (45|. We also show the bias in the parameters M and v when achieving the 
effectualness. For each pair, the top number is Initial LIGO while the bottom number is Advanced LIGO. The sign of the bias 
is such that in all cases the fiducial EOB templates slightly overestimate the total mass M and underestimate the mass ratio 
v of the improved EOB signal. 





Effectualness AM/M Av/v 


Faithfulness 


(1.4, IA)Mg) 


0.999 0.98% -1.63% 
0.999 0.98% -1.63% 


0.992 
0.995 


(1.38, 1.42)M Q 


0.999 0.96% -1.60% 
0.999 0.89% -1.49% 


0.992 
0.995 


(5,5)M 


0.997 1.32% -2.12% 
0.999 2.06% -3.47% 


0.973 
0.976 


(4.8, 5.2)M 


0.999 2.42% -4.08% 
0.999 2.11% -3.54% 


0.973 
0.976 


(10, 1O)M 


0.999 2.70% -4.62% 
0.999 2.59% -4.39% 


0.974 
0.962 


(9.5, 10.5)M Q 


0.998 1.40% -1.94% 
0.997 2.67% -4.54% 


0.974 
0.964 


(15,15)M Q 


0.995 4.80% -9.98% 
0.999 2.49% -4.23% 


0.987 
0.973 


(25,25)M 


0.995 4.95% -12.6% 
0.994 3.00% -5.56% 


0.982 
0.985 



frequency is artificially extended to a higher frequency. However, this might require a tweaking of the upper frequency 
cutoff depending on the noise spectral density of the detector [36| and the mass ratio of the system, and either the 
extension to unphysical values of v [35|, H(| or the inclusion of a p4PN term in the template phase [35| calibrated to 
the numerical simulations. We believe that a better alternative for heavier systems are the EOB templates calibrated 
to numerical relativity simulations [H, H3, HI, S S3, 0, S3, S3, H, HBJ . The most recent EOB models are in near 
perfect agreement with the most accurate numerical simulations to date, although only a small number of systems 
corresponding to different mass ratios have been studied so far. Nevertheless, a physical model with physically 
meaningful parameters is a far safer bet as search templates unless, of course, if the model in question is not in 
agreement with the waveform predicted by numerical relativity. So far, the EOB is the best physical model we have 
and this is what we recommend be used to search for binaries with mass es g reater than about 12Mq. 

In this paper we adopted the preliminary, fiducial EOB model of Ref . [39[ > because it is the EOB model currently 
available in LAL and it is used for searches by Initial LIGO. For completeness, here we quantify the closeness between 
the EOB model used in this paper and a most recent improved version of the EOB model [45| (which is similar to 
the one of Ref. [Hj]). The latter was calibrated to longer and more accurate numerical waveforms generated by the 
Caltech/Cornell pseudo-spectral code [8(|. Reference [45| found that the faithfulness of the improved EOB model to 
these highly accurate numerical waveforms is better than 0.999. In Table HVl we show both the faithfulness and the 
effectualness of the EOB model [39| to the improved EOB model [45| using noise spectral densities of Initial LIGO, as 
well as the bias in the parameters M and v when achieving the effectualness. The search for effectualness in this test is 
done continuously in the parameter space, instead of using a template bank. Although there is some systematic trend 
in the numbers due to the difference in the EOB models, the main result is that the faithfulness and the effectualness 
are always better than 0.97 and 0.995, respectively. Assuming the numerical waveforms of Ref. [45[ are exact, the 
EOB model of Ref. [3^ used in this paper is accurate for detection purpose with a loss of event rates of ~ 10%, and 
may cause ~ 10% bias in estimating the mass parameters. 

In this study we considered PN waveforms in the so-called restricted PN approximation. Restricted waveforms 
contain only the second harmonic of the orbital frequency. Inclusion of other harmonics is necessary, especially 
when a binary is arbitrarily oriented with respect to a detector and the component masses are dissimilar. Recent 
studies [60l . IrJlT ] have shown the tremendous advantage of including these other harmonics in the GW templates. 
Therefore, it is necessary that a future effort undertakes a study similar to this, but includes all the amplitude 
corrections. Furthermore, Ref. [46| has shown that by supplementing the PN results by the available test particle 
results up to 5.5PN improves the match between the EOB models and numerical relativity simulations. This can be 
expected to lead to further improvements in the results obtained here in the future. 
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